Green's function technique for a two-electrode mesoscopic system under bias 
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We present a Green's function technique for studying the nonlinear conductance of a nanocontact 
system with two electrodes at different chemical potentials. The retarded Green's function for a 
single-impurity Anderson model with two reservoirs is obtained in terms of a 5 x 5 matrix in which the 
effect of bias is contained. A complete set of basis vectors for the single-impurity Anderson model 
has been provided before formulating the Green's function. Finally, we present a self-consistent 
method to fix the undetermined quantities existing in the matrix elements for the retarded Green's 
function. 

PACS numbers: 72.10.Bg, 72.15.Qm, 85.75.Mm, 73.23.-b, 73.63.-b, 75.76.-f-j 



I. INTRODUCTION 

The single-particle Green's function is a basic tool for 
studying the correlation effects in many-body systems. 
However, calculation of this function in the case of a 
strongly correlated system is not usually successful be- 
cause such a calculation generally requires nonperturba- 
tive treatment. Moreover, the possibility of a successful 
calculation becomes even less when the system is under 
steady-state nonequilibrium (SSN) conditions. Recently, 
some mesoscopic systems having nanocontact features 
have attracted considerable interest. Such systems in- 
clude a quantum point contact llHj|, a single electron 
transistor with a quantum dot [Sf-ls ll , and magnetic atom 
adsorbed on a metallic surface [9h12!|. Experimentalists 
measure the dl /dV in these systems for a given V; I and 
V denote the current and the bias, respectively. Because 
of the bias, resonant tunneling occurs unidirectionally 
in the ground state. This unidirectional tunneling causes 
the system to be out of the linear response regime. There- 
fore, nonlinear conductances observed for the abovemen- 
tioned systems may not be explained by a linear response 
theory. The purpose of this study is to provide a Green's 
function technique that is suitable for studying the non- 
linear conductance of a mesoscopic system operating un- 
der bias. 

A well-known formula for a steady current passing 
through a small interacting region that is connected to 
charge reservoirs having different chemical potentials was 
proposed more than a decade ago [iSJ, 11J| . This formula 
is given by 
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for proportional lead functions, i.e., T^(u}) ex T^{lij) — 
27r^^/9^(aj)[Vj^P, where a indicates spin, L and R de- 
note the left and right reservoirs, respectively, p^{u!) is 
the density of states of the right metallic reservoir with- 
out inter-electron interaction, and Vj^ is the strength of 
hybridization between the state d of interacting region 
and the state k of the right reservoir. Further, /l(w) in 



Eq. ([T]) denotes the Fermi distribution function of the 
left reservoir. Since unidirectional tunneling at ground 
state is a unique characteristic of the SSN, the effect of 
bias must be appeared explicitly in the retarded Green's 
function in Eq. ([IJ. However, to the best of our knowl- 
edge the steady-state Green's function for the quantum 
impurity system has not been obtained yet. 

The purpose of this study is to present a method for 
determining the steady-state retarded Green's function 
appearing in Eq. ([T]). To realize this method, we adopt 
the retarded Green's function expressed in resolvent op- 
erator form in the Heisenberg picture [15| : this function 
is given by 



^G'+,(z)-(q,|(zI + zL)-1|c,,), 



(2) 



where z = —iuj + rj and Cia indicates annihilation of a 
fermion with spin cr at a state i. The symbols I and 
L in Eq. ^ denote the identity operator and the Li- 
ouville operator, respectively. The latter is defined by 
LO = [n, 6] ^ no - on, where n and 6 are the 
Hamiltonian operator and an arbitrary operator, respec- 
tively. Although we are familiar with the resolvent form 
in the Schrodinger picture, 



G+H = (0.|(c.-|-ir/-H)-i|(/),), 



(3) 



where (j)j is an eigenstate of n, we employ Eq. ^ in this 
study because a dynamical approach using operators is 
more appropriate for describing electron hopping under 
steady-state conditions. 

This paper is composed as follows: In section H, we 
present a systematic method to determine a complete 
set of basis vectors for calculating the retarded Green's 
function given in Eq. ^. A single-impurity Anderson 
model with one- and two-reservoir is studied as an exam- 
ple. In section HI, we express the retarded Green's func- 
tion for a two-reservoir system under bias in terms of a 
finite-dimensional matrix and discuss the physical mean- 
ings of the matrix elements. We present a self-consistent 
method to determine unknown quantities in section IV 
and finally give conclusions in section V. 



II. DETERMINING COMPLETE SET OF BASIS 
VECTORS 

The first step in the calculation of the retarded Green's 
function in resolvent form, Eq. ^, is ensuring a com- 
plete set of basis vectors. It should be noted that the 
resolvent form does not provide information about the 
orthogonal basis vectors. Therefore, the Lanczos algo- 
rithm [ig, |l7[ in the Schodinger picture and the projec- 
tion operator technique [18|, [i9| in the Heisenberg picture 
are used to obtain static and dynamic orthogonal ba- 
sis vectors, respectively. However, employing the Gram- 
Schmidt orthogonalization procedure to obtain orthog- 
onal basis vectors for a nontrivial system causes signifi- 
cant complications. Hence, we propose a new and simpler 
technique for obtaining the basis vectors. 

The retarded Green's function in Eq. ([2]) can be ob- 
tained by calculating the matrix inverse (M~^)ii, where 
M is composed of the elements 

Mij = z6ij - {iLBjlei) = z5ij + (ejlihe^), (4) 

once a complete basis set, {e£\£ — 1, • • • ,oo}, spanning 
the Liouville space is given. We select ei — Cda to ob- 
tain G'^citri^)- "^^^ inner product in Eq. ^ is defined by 
{ek\ei) = ({efe,e|}), where the angular and curly brack- 
ets denote the statistical average and anticommutator, 
respectively. We now present a systematic method for 
determining the basis vectors. For this purpose, express- 
ing the Green's function operator in a diagonal matrix 
as 
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is the most crucial step. Once this form is achieved, the 
Liouville space of G^^(z) will be spanned completely 

by the linearly independent components of vectors \A) 
and \]3). To obtain the form in Eq. (O, we expand 
the Green's function operator G = [z\ + «(L/ -I- Lc)]^^, 
where L/ and Lc respectively represent the Liouville 
operators for the isolated part "Hi and the connect- 
ing part T-Lc of the Hamiltonian; The operator is ex- 
panded in powers of Lc using the operator identity 
{A + B)-^ == A-^-A-^B{A + B)-^, where A = zl + iLi 
and B ~ ihc- Then, after resumming the expansion, the 
retarded Green's function can be expressed as follows: 



iGJ^^[z) = {cia\Gi\cj„) - {ci„\Gi\iLcGiCjc 
+ {ciaGiiLc\G\iLcGiCj„), 



(6) 



where Gi = [zl + iL/) ^ . Equation ^ can be written 
in the matrix form as 
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where G 
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\iLcGiCja), and 



Gi -Gi 
the superscript T denotes the transpose. Using the lin- 



ear transformation U 



I -Gj 





^^ I, where Glt 



Gi/[G — Gi], one can diagonalize G. Equation ([7]) is 
then modified as 

^GtM = {{c..A ($.|)Gd(|c,.) |<i>,)f , (8) 



where Gn = ( ^ ) and 

' G 



|c,,) + (G,-G)-iG,|$,) 



(9) 



Equation ([8]) is the desired form and the linearly inde- 
pendent components of vectors \c.j^) and |$j) completely 
span the Liouville space of Gj^j^{uj). In conclusion, the 
systematic method for collecting the basis vectors in- 
volves determining all linearly independent components 
comprising the vector \cja) because |$j) is contained in 
this vector. 

Wc will subsequently demonstrate the determination 
of the basis vectors by means of an example. We are in- 
terested in a mesoscopic system with a mediating Kondo 
atom between two metallic reservoirs. This system can 
be described by a single-impurity Anderson model with 
two metallic reservoirs, whose Hamiltonian is expressed 
as 
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where e^, e^, Vkd, and U indicate the energies of an elec- 
tron of momentum k in lead, the level of a mediating 
atom, the hybridization between the atom and lead, and 
the on-site Coulomb repulsion at the atom, respectively. 
Because the two-reservoir Anderson model is a straight- 
forward extension of the one-reservoir model, we consider 
a one-reservoir Anderson model as an example. Its iso- 
lated and coupled parts are given by 



ni = 



] f-kcl^Cka + E '^dcl^Cdcr + Uud^Udl, (10) 
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and 
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We first demonstrate the calculation of ^d — *Lc(I + 
iLi)~^Cd'\ and then of \cda) for T-ij and He above. It 
can be clearly observed that (I + iLi)^^Cd-x yields only 
two operators, Cdi; and ridxCdi;- Applying Lc to these 
operators changes the index d{k) into k[d) in fermion 
operators and Udi into j^, . Therefore, the components 
of <I>d are given by 



$ 



d^{ck\, ndiCkf, Jd^Cdf), 



(12) 



where fc = 1, 2, ■ • • , oo and 



■Idi 
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We now focus on the vector \cdcr) in Eq. ([9]). The oper- 
ator L$d = [H,ck{\ + [H,ndiCk{\ + [H,JdiCd-\\ yields Cd-\ 
and Cfcf from the first commutator; nd^Cdt, JT-diCfc-t-, and 
Jd\Cki; from the second commutator; and Jd\Cd-\, Jdi'^k^^ 
Jdi'^dt, Jdi'^diCdf, and jj4,?^dtCdt from the third commu- 
tator, where 



^di 



J2(^kdclcki + V,:dclcdi). (14) 



The last operator jji ^^dfCdt is a vanishing one, and the 
operator Jd\''T'diCd-\ is dynamically equivalent to JdiCd-f- 
Thus, the linearly independent components of (I + L)$d 
are classified into two groups, one involving Ck^ and the 
other involving Cd^. These components are 

(I + L)*d = (cfct, n-diCfct, Jd^Ckt), (15) 

for fc = 1, 2, • • • , oo, and 

(I + L)$d = (cdt, ndiCdt, jJ^Cdt, J^'^Cdt). (16) 

Finally, we apply L^ , which is equivalent to the re- 
peated application of Lc, to the operators in Eqs. (fTSJ) 
and ([TBI) . It can be clearly observed that the multiple 
application of Lp to Ckf, Cd^, and Udf simply repro- 
duces the existing operators and Jdi'^kt- The remaining 
linearly independent operators arise from L^jJ^, where 
n = l,2,---,oo. Thus, one may classify the basis vectors 
into two groups. 



Set of Cd^ 



{cdt, ndiCdf, Jdi^dt, i^Cdt, 
i^C&dt, (LcJ7i"d4)cdt} (17) 



and 



Set of Ckf 



{cfct, ndiCkf, 
ih-jl)ck^}. 



JdiCkt, 



J^Cfct, 



(18) 



The operators shown in Eqs. ([TT]) and ^TE\\ describe all 
the possible linearly independent ways of annihilating an 
up-spin at site d in time t. In other words, these operators 
completely span the Liouville space of Cdf{t). 

It is important to identify the meaning of the basis 
vectors (Lg,j^)cfc^ and (Lg.jJ^)cd-|- for later discussion. 

Application of Lc to jj, gives rise to 



'^CJdi 



1 k 



2i^|VkdP4^Cd; 



and 

Lcj/'i = XI Zl VldVCdcl^Ck^ -J2Y1 ^W^kd^^Cu. 
Ik Ik 

These operators represent a round trip of a down-spin 
electron between the mediating atom and electron reser- 
voir. Therefore, (L^j^,) depicts n-time trip of a down- 
spin electron between the mediating atom and the reser- 
voir without coming into contact with an up-spin electron 
at the mediating atom. This phenomenon rarely occurs 
in reality, and we neglect these basis vectors. 



III. RETARDED GREEN'S FUNCTION FOR 
TWO-ELECTRODE SYSTEM 

We identify the basis vectors that do not play an im- 
portant role in describing the dynamics Cd^{t) in the 
Kondo regime. We first consider the basis vector nd^Cd-f 
in Eq. (ITT)) . Since n'^^ = Udi, this basis vector describes 
the dynamical processes involving all higher orders of 
Coulomb interaction U. Therefore, in the Kondo regime, 
ndiCdf must be eliminated. The operator j^/ndiCd^ in 
Eq. (fT7|) also represents the process that costs energy 
U or more. Therefore, we eliminate the basis vectors 
{'LQJZndi)cd^ in studying the Kondo regime. We further 
eliminate the basis vectors (L^jT )crf-|- because multiple 
round trip of an up-spin electron without coming into 
contact with a down-spin electron at the mediating atom 
are practically rare. These approximations are valid for 
describing the Kondo processes. Now, we discuss on the 
basis vectors combined with Ck shown in Eq. p8)) . These 
basis vectors contribute to constructing the self-energy. 
It is sufficient to choose Ck and nd^Ck^ in Eq. p8|l to con- 
struct the self-energy. The basis vectors (L^j^)cfc^ are 
eliminated by the same reason eliminating (L^j^,)cd-x i^ 
Eq. (|17p and the contribution by jji c^^ is much smaller 
than that by n^^Cfc-i-. 

Finally, we obtain the reduce Liouville space con- 
structed by the basis vectors 

i^df, Jdi^dt, itiCdt) and (cfet, ndiCkf), 

where k = 1, 2, ■ • • , oo for describing the dynamics of the 
single-impurity Anderson model in the Kondo regime. 

Since the resolvent form of Gj^^(z) is given by the 
Laplace transform of the coefficient ai{t) in the expan- 
sion Cd^{t) = Xli^i 0'dt)^ii where ei = c^^, the equality 
ai(i) = G'^^Jt) must be obtained. This equality is valid 
when all the basis vectors except Cd-[ are orthogonal to 
Crf-f-. This orthogonality condition is satisfied by introduc- 
ing the expression 5A = A — {A) for the operators n^t 
and j£. Thus, we finally obtain the reduced Liouville 
space spanned by the basis vectors 

(cdt, Sj'^^Cd^, ^Jdi^di:) and (c^^, Sud^Ck-x), 

where k — 1,2,--- ,oo, for the one-reservoir Anderson 
model in the Kondo regime. In practical calculations, 



we use the normalized forms of the basis vectors. The 
basis vectors given above are the same as those used in 
the previous study [20|, in which the basis vectors were 
obtained intuitively. The Liouville space for the two- 
reservoir Anderson model is a straightforward extension 
of that for the single-reservoir model. The basis vectors 
are given by 



(c" 



fcf 



SudiCf^-f), where k = 1,2,- ■ ■ , cxo; and 



-■R, 



i^Jdi Cdt, SJdi Cdt, Crft, Sj^^ Cdt, Sj^^ Crft) and 



{Sudic 



^kf) 



where k — 1,2,- 



,00, 



where the superscripts L and R denote the left and right 
metallic leads, respectively. The matrix M of Eq. (|4]) is 
expressed as 



Mll MdL 
M = I MLd Md Mm 
Mdn. Mr„, 



(19) 



where Md is the 5x5 block that is constructed by the 
basis vectors 



(St' 



'di Cdt, 



'^Jd/'Cdt, 



Cdt' f^Jdf^^dt 



Sjlfodt); 



Mll, the 00 x 00 block constructed by the basis vectors 
{cj^f, Sudicj^-t-), where fc = l,2,---,oo; and M^l, for ex- 
ample, a 5 X 00 block constructed by the basis vectors 
given above. 

The retarded Green's function G^^^{uj) is obtained 
by calculating the matrix inverse (M~^)dd, where the 
subscript dd indicates the element obtained from Cd^ 
in both row and column. Calculating the inverse 
of an infinite-dimensional matrix M is a nontrivial 
problem. For this purpose, we perform matrix re- 
duction by using Lowdin's partitioning technique [2ll . 
|22| . We consider an eigenvalue equation MC = 
0, in which the column vector C — (C^ C^ C/j)-^, 
where Cl, C^, and Cr correspond to the row vectors 
(Cfe^, Sndi4^), (SJdi^Cd^, SJdi^Cd^: Cdt, ^J7i"cdt, Sjjfcd^), 
and ((5nd4.c£,,cS.), respectively, and is a column vec- 
tor with zero elements. After eliminating C^ and 
Cr, we obtain an equation (M^ — MLd^LL^dL ~ 
MRdM-^]^MdB.)Cd = MrCd = 0, which leads to 

M, ^Md- MLdMllMdL - MRdM^]^MdR. (20) 

The last two terms in Eq. (PP)) contribute to the self- 
energy of the retarded Green's function. We finally ob- 
tain Mr for the two-reservoir Anderson model as follows: 
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Here, w' = u] — ed — U{nd\), where e^, U, and {nd\,) denote 
the energy level of the Kondo impurity, Goulomb interac- 
tion, and the average number of down-spin electrons oc- 
cupying the level e^, respectively. The derivations of ma- 
trix elements are given in Appendix A. All the matrix ele- 
ments, except the eight [/-elements, have additional self- 
energy terms iT^mn = /3™„(iE^(tj) -I- iT,§{uj)) = 2/3„„A, 
where I]^(a;) = J^k I^mIV(^ - ^k + iv) = -*A^ for a 
fiat wide-band and A ee (A^ + A^)/2. A is used as 
a unit of energy. The lead function r^(a;) in Eq. (1) 
is equal to 2A^ in this study. The coefficients Pmn are 

/3ll = IC'^P, (322 = IC'^P, /333 = 1, /344 = krP, /355 = 
\Q"f, /3l2(21) = [C^*C'^](*\ /3l4(41) = [C^*C''](*\ 
/3l5(51) = [C'^*^'']^*^ /?24(42) - [C'^*^^]'*), /325(52) = 

[e+^*C^]/*\ and /345(54) = Kr*C'']^*^ Where ^J\ for 
example, is expressed as 



cTi _ 
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We show the explicit expressions of pmn in Appendix B. 
The matrix elements represented by 7 are written as 



-L(fl) ■+L(Rh 
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We show the detailed derivations in Appendix A. The 
retarded Green's function iC^^Auj) — (M~^)33 may 
explain all the features of nonlinear conductance ob- 
served in a quantum point contact [l|-[j], a single electron 
transistor j^-Q, and scanning tunneling microscopv on a 
magnetic atom adsorbed on a metallic substrate |9l-[l3|. 

In Fig. 1, we present the graphical illustrations of 7 
on the basis of the corresponding operator expressions, 
which are third-order processes of hybridization. Figure 1 
depicts the Kondo processes in which a down-spin indi- 
cated by first enters the Kondo impurity and forms a 
singlet with an un-spin electron indicated by (2). Then, 
it performs an exchange process or singlet hopping. Fig- 
ure 1 (a) describes the spin movements in ^ll (left) and 
iRR (right) in equilibrium, and Fig. 1 (b), those in ^lr 
{+) and 7j (— ). Each 7 is composed of two terms, i.e., 
spin exchange and singlet hopping. The latter is a unique 
property of a two-reservoir system. It is noteworthy that 



FIG. 1; Motions of electrons indicated in 7. (a) Kondo pro- 
cesses described by 'jll (left) and '^rr (right). The solid dot 
indicates a Kondo impurity. The up- and down-spin motions 
denoted by (2) on the same side represent the exchange pro- 
cess while the up- and down-spin motions denoted by @ on 
different sides indicate singlet hopping, (b) Kondo processes 
of connecting mechanism: '^lr (+ sign) and 7j (— sign), (c) 
Unidirectional resonant tunneling of a singlet establishing cur- 
rent flow. 



7ll and ^rr represent the degree of Kondo coupling, 
whereas ^lr and 7j represent the resonant tunneling be- 
tween two reservoirs. 

The effect of bias is contained in the Fermi distribution 
function and is observed in the calculation of the expec- 
tation values of 7. In the ground state, resonant tunnel- 
ing from the reservoir at lower chemical potential to the 
one at higher chemical potential is prohibited. There- 
fore, the second part of Fig. 1 (b) vanishes. The motion 
that establishes current flow in both ^lr and 77 is uni- 
directional, as shown in Fig. 1 (c), which leads to the 
condition ^lr = IJ- However, when the second part of 
Fig. 1 (b) does not vanish for certain reason, the condi- 
tion 7ifl > 7j is obtained. In both conditions of 7, even 
the state of very low bias is quite different from that of 
zero bias at which 7,7 = 0. 



IV. SELF-CONSISTENT METHOD 

The values of 7 and the fluctuations in the denom- 
inators of Eqs. (P^ — ([^S)) may not be determined by 
direct calculation. These values can be determined by 
employing^ the self-consistent method, which was used by 
Nagaoka[23| for studying the conventional Kondo prob- 
lem. However, we do not try to determine these values in 
this study. Nevertheless, in the following study, we con- 
sider them as free parameters in order to explain the ex- 
perimental observations for a quantum point contact [l|-y| 
and scanning tunneling microscopy [2-[l3|. 

Finally, we suggest a self-consistent scheme that is par- 
ticulary useful when the application of a magnetic field 
disturbs the particle-hole symmetry. Because the ma- 
trix elements U j^ and coefficients Pmn contain undeter- 
mined quantities {uda}, {Jdc ' }: and (jj^ ' ), and these 



quantities are also expressed by the retarded Green's 
function such as 



{Uda} ^ -- j 



r^(u;) + r«(a;) 
X ImG+,(c.)dc.; (26) 



07/) = - 



[!L{uj)~fR{uj)]T{u) 



1 r dco 

X imG+,,H = -07/); 



(27) 



.+i(fl)\ 



duj 



Um = I ^/LwHr^(^n^)ReG+,,(a.), (28) 

where f(a;) = r^(w)r^(a;)/[r-^(w) -H r«(w)], one can 
construct a self-consistent loop to obtain the retarded 
Green's function as follows: 



.(")i 



(n.,)(°),(j,f^«)(°) ^ G+V°)(a.) ^ p-"iu) 

{n,iy'\ 0.f ^^y^^ ^ Gtll\^) ^ P?\^) ^ • • • 

The expressions in Eqs. (|26l) and (P7| are valid only 
when r^(LL') oc r^(a;). We report the results for the 
magnetic-field-induced peak splitting of the differential 
conductance in a single-electron transistor with a quan- 
tum dotp-i8] in terms of this self-consistent scheme in a 
separate study. 

V. CONCLUSIONS 

In conclusion, we have presented a systematic method- 
ology for determining basis vectors spanning the Liouville 
space, which is the most crucial step in calculating the 
retarded Green's function using the resolvent operator; 
further, we have suggested a procedure for calculating 
the retarded Green's function. The operator method pre- 
sented in this study has several advantages: (i) A com- 
plete set of basis vectors cab be determined systemati- 
cally, and (ii) the physical meanings of the basis vectors 
represented by the operators are apparent, allowing the 
identification and removal of the unimportant basis vec- 
tors for a particular parameter regime and subsequent 
construction of a reduced Liouville space. The on-site 
retarded Green's function for the single-impurity Ander- 
son model with two metallic reservoirs is expressed as an 
inverse of a 5 x 5 matrix. We discussed the characteristics 
of the matrix elements under the application of a bias. 
Finally, we suggest a self-consistent calculation method 
that is useful when the system under consideration does 
not show particle-hole symmetry. This operator formula- 
tion for the retarded Green's function may be appropriate 
for studying the dynamics of a many-body system under 
steady-state nonequilibrium. In a separate study, we ob- 
tain the spectral function and differential conductance 
for a specific system that has a Kondo impurity between 
two metallic reservoirs under bias. 
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Appendix A: Calculation of matrix elements 

Galculation of the matrix elements of Eq. (19) for the 
single-impurity Anderson model with one reservoir in 
which ^LR and 7j do not appear is shown below. One 
can obtain the expressions of 'yLR and 7]^ via the same 
manner shown here. The block M^^i^^^) comprises two 
diagonal blocks whose elements are the same and given 
by —iLo + iei + 0"'", —iuj + ie2 + 0^, and so on. Since the 
calculation is very simple, we skip it. 
(1) Matrix elements of the block M^^: 
Nontrivial elements of M^d are 

-{{iL{cki:6ndi),6jf^cli^}), i.e., ({i[H,CktSndi],jf^cli^}) 
= {{i[H, Ck^]Sndi + Ck-\i[H, 5ndi\,5j^^c\^}) 
= -iek{{ck^5ndu5i^^c\^}) - iVkd{{cd^5ndu5i^^c\^}) + 

({cfet'5j7i''^i<ii4t})- '^^^ ^'"'^^ *^^"^ 

{-iek){{ck\5ndi,5jj^c\^}) = ^ (-Jefc)([ndi, jljcfct^^) 
should vanish because the up-spin dynamics of this form 
is not allowed. The third term also vanishes because 

({cfetJli: '^jli^tD = 0Ti'^jli{cfct>4t}) = 0- The second 
term, however, must be calculated rigorously, since it 
is a hybridization term that is related to the Kondo 
process. The operator {cdi;5ndi^T5j^,(y^A is expanded as 

Cd^[5ndi,5j^^c\^] + [5j'f^c\^,Cdi;]5ndi + 2cd^5j^^cl^5ndi. 
Each term is rewritten as Cd-\[5nd\,T5j^.c}^^} ~ 
CdtKi,.?l^]4t == Cdt^tKi, jJJ = (l-nrft)[nd;,jJJ and 
[^Jli^^dv (^di]5ndi = 5i^^[c\^, Cdi]5ndi = 5i^^{l - 2cdi;c\^) 
5nd\,. The last expression —2cd\5j^^Cj^^5ndi cancels the 
third term above. Hence, — ({iL(cfcf(5nd4,), jJ c]j^}) — 
«Vfed((l - nd^)[ndi,5j^^] + SjJ^Sudi) = 

iVkd{{^~ndt)[ndi,j]:^] + 5Ki,iJJ + Sj];^Sndi) = 
i^iil 2nd,)[nd,,jl] + {Snd,,6jl}) = ^-^{{{1 

2ndt)K;, jJJ) + (1 - 2(nd^))(jT )}. 
Therefore, the matrix elements 



ied{{cdt, 4t}) + ^U{{cdtndi, cl^}} 
= led + iU{{cd^ndi, c\^}) = ied + iU{{cdt, cli^}ndi) 
= ied + iU{ndi), and 
(B):-{{t[H,CdtSjl]AcdtSjlV}) 

= -mH,Cdt]Sjl,Sjlcl^}) - {{cdAH,Sjl],5jlcl^}). 
The first term of (B) is rewritten as 
{{i[H,Cdt]SjJ^^,SjJ^^cl^}) = ({(-iEk^kVkt -*edCdt - 
iUcd^ndi)SjJ^, {cd^Sjf^y}) = -ied{{cdtSjJ^,SjJj^cl^}) - 
iU{{cdtndiSjJ^,SjJ^cl^}). Applying the decoupling ap- 
proximation, ndiSj]^, — {ndi)Sjj^i, to the L'^-term above 
gives rise to a form of squared norm for the first term of 
(B), i.e., {{i[H,Cdt]Sj]:^,Sj]:^cli^}) = [-ied - iU{ndi)] x 



\cdt<>j, 



dV 



The second term of (B), however, cannot be 



written in the form of a squared norm, because [H, j^^A 

j^^. Therefore, we neglect it. 

(2-2) Matrix elements Uj^ : 

These are given by the inner products 



oc 



-{{iLcdt,SjJ^cl^}), i.e., -{{i[H,Cd^],SjJ^_cl^}) 



tJ 



= ^'^d{SjJ:^){{cd^, cl^}) + iU{{cdtndi, SjJ^cl^}) 

+ ^Ekl4*d(^jli)({ckt,4t}) = ^U{{cd^nd:„Sjlcl^}). 

However, {{cdtSudi, SjJ^cli^}) = {cdt[6ndi, 6jJ^cli^]) + 

ii^Jdi^d^'^d^]^ndi) + {2cdfdjJ^cl^^Sndi), where 

the operator in the first term is rewritten as 

Cdi:[Sndi,SjJ^cl^^] = Cd^indi, jj{\cl^ = Cdt4t["di'.?JJ = 
(1 — nd^)[ndi, jji], while the second term is given 
by {[SJdAv('dt]Sndi) = {Sj^^^[cl^,Cdt]Sndi) = 
{5jd^6ndi) - (2(5jJj_Cdtcj,^(5ndi). The last one can- 
cels the third term of {{cdfSudi, ^Jdi'^d-f}) above. Hence, 
{{cdi:Sndi,Sjf^cli^}) = ((1 - ndi:)[ndi,jj^] + Sjf^Sudi) = 
{{l-ndt)[ndi,j^l] + U^ndi,SjJ^] + SjJ^Sndi) = 
i((l - 2ndt)[ndi,J^^] + {Sudi^Sj^^}) = i{((l - 
2ndi:)[ndi,jf^]) + (1 - 2{ndi)){jJ^)}. Therefore, the 
matrix elements denoted by C/jt is given by 



U 



JT 



-{{^^Cd,,Sjlcl,}){i5jir)-'^^ 



iU, 



= —[{{l-2ndt)[ndi,jf^]) 

+ (l-2(n,,))0-J^)]((<5jJ,)2)-V2. 



(2-3) Matrix elements Jll{rfi.)- 



-{{iL{ck^Sndi),6jJ^cl^})/[\\ck-f5ndi\\ x ||Qt(5jJJ|] = The inner products -({iL(cdt(5jJ^), (5.7^^4^}) gives rise 



Vkd^'J, where ^J = Hi^i^ - 2ndf)[ndi, jj^]) + i{l - 
■2{ndi)){3l)}{{Snd,r)-'/H{Sjir)-'/\ 

(2) Matrix elements of the block M^: 

(2-1) Diagonal elements: 

By using the commutator expression i[H,Cd-\-] = 

~i Ek ^kd^kt - iedCdt - iUcd-^ndi-, the diagonal elements 

of the block M^; are given by 

(A): -{{i[H,Cd^lc%}) - ^Ek^kd({ckt,4^}) + 



to the expressions for Jll(rr)- Since 

({*L(c,t<5j,\),<5j±4^}) - {{^[H,Cdtm^,Sj^^c\^}) + 

{{cd-fi[H,j]]],Sjj^:cliA), the first term can be expanded 

as {{i[H,Cdt]Sj];^,Sj^^cl^}) = -ied{{cdtSjJ^,Sj^^cl^}) - 

iU{{cdfndiSjJ^, Sjf^^}) - i Ek KdiicktSjJ^, ^JdAt^^ 
The first and second terms of this expression should van- 
ish because they are not written in the form of a squared 
norm. Since the third term describes the Kondo process, 
we treat it rigorously and have the following expression: 



= -*Ek^kd(bdl.it]'='=tCdT)- The second term of 
({iL(cdt^jl;),5j^;4^}), i.e., ({cdti[7J,jJJ,5jj=^ct^}), is 
rather comphcate because it contains the commutator 
i[i/, jS]. This commutator is expanded as i\ii^f^A = 

^^d3% T («) Efe efe(VfcdC^Cd; ± Vkdc\^CkC) ± t^^dfJ^; ± 

(*)EM(VfcdV£d4^Q4, ± yudYuAiCn). where (i) is ap- 
phed only to the lower signs and the prime in sum de- 
notes k ^ L The first two terms on the right side are 
cancelled if we assume Ck ~ fd, and the [/-term is ne- 
glected in the Kondo regime. The last term, second or- 
der hybridization, describes the round trip of a down- 
spin electron between the Kondo impurity and the lead, 
which is equivalent to (j^)^. Since we consider only 
single trip of a down-spin electron in this study, we ne- 
glect the term involving [H,j^A. Thus, the matrix el- 

ement (Md)2i/i2 = — , ^, ]: ^ = ±7, where 



7 = ± 






for the single-reservoir An- 



derson model. ^LL(RR) for the two-reservoir Anderson 
model is given by 



ILL(RR) 



iY" V* (c^ r^ -^c^c^ Mr^^'^'^ j+^'-^'h) 

((^J-^(«^)2)l/2((5_,-+^(«))2)l/2 ■ 



Appendix B: Expressions of /3m„ 

We obtain the expressions of /3m„ in terms of the def- 
initions and the expressions of S,^ given in Eq. (1221) . 
The coefficients jSmn in front of the self-energy function, 
i.e., iY,mn{i-^) = Pmn{i^o{^) + «S^(w)), are symmetric 
in exchanging their indices. By using the operator iden- 
tity [jT-i^, jj] = T^Jiii we write a different expression for 
Ke[U^t\ i.e., Re[U^t"] = {U/2t), where 



1/r 



[iO-"^"") T 20;f • Vt)]/ ^/((^jjp)^, 



which will be used in the expressions of Re[/3mTi]. The 
final forms of Re[/3m„] are given by using ((^4)^) « 1/4 
and T = 2 that is the value at the atomic limit. Since 
(j+'^(''^) < and (j-^) = -Oa""), an inequality, 

R-e[/3i2(i4)] < Rc[/3i5] < Re[/3ii(55)] < Re[/325(45)] < 
Re[/322(44)], exists. We obtain the expressions of Pmn as 
follows: 

Real parts of /?,„„: 

Re[/3n] = MCl*Cl] = 



L\2 



((jr)-2Wjr))'+(i-2('^^i))'(j4 ) 









4{(<5nu)"> 



1 (1-2(»,^))^(j7/)^ 

Only the index L changes to R for Re[/355]. 



Rc[/?22] = RcKre^] = 

(07/) - 2KtJ7/))2 + (1 - 2{n,ClY{jt,^) 



L\2 






1 {J~2{rH^^f{g^ 

^ {('5ji'') = > 

1 (l-2(«,4))"(j+7)" ' 

4 ^ /i'Ai-+-^12 



4((5«.4.)2) 



Here, again, the index L changes to R for Re[/?44] 
Re[/3i2] = Re[eZ*e^] = Re[/32i] - 

+ (1 An,,)) 0,, )0,^ )J ^___^__^ 

1 ^ (i-2K^))^077)0-+7) 

Rc[/3i4] = Rc[C*eA] = MPii] = 
'iijt') 2(n,,j+^))(-07^) + 2(n,,j7^)) 
+ (l-2(n,,))207/)0-+«)" 



^ii^lIi^^iiSjt^ 



+ 



(1-2(«,4,))^(j7)0-+7) 

(i-2(«.^))"o-7)o-+") 



4{(5n.4)2) 



Re[/3i5] - ReKl*^;,] = Rc[/35i] = 
"(O'i") " 2(n,,j+/))(0-+/) - 2(n,tJ+/)) 

+(l-2(n,,))207^)07^)l ^M^-^niii. 



yu^j^n/aV^ 



(i-2(».4.))^(j-7)o--7) 



/((^j^fH/{(V¥> 



1 , (l-2{n.^))^0-7)077) 



4((5«.J-) 



-K 



Re[/324] = Re^reA] = Re[/342] = 
(2(^.tJ4^) - (J4''»(2(n.tj4^) - (j^«)) 
+(l-2(n,7)20+^)0+^) 



[4((g«u)')]' 



a/u^jS^v^o^O 



-h 



(i-2(,».^))^(j+7)0+7) 

(l-2(n.^))^0+7)(j+«) 



4{(fc>4)^) 



Re[/325] = ReKrCj,] - M^^2] 



{-{j^) + 2(n,,j^^))((j+«) - 2(n,,j+«)) 



+(l-2(n,,))2(,+^)(,-«) 



[4{(Jn.J^)]-i 



1 I (i-2(».4,))^0-+/)(j7/) 

1 , (l-2(nu))^0+x^)0u") 



4((fc.4,)^) 



Re[/345] = RcKA^Ck] = Rc[A54] = 



^((^jiv>v^((^0^ 



4{(5n.^)^> 



1 I (i-2(».,))^(j + ")07/) 



Imaginary parts of /?„„: 



(l-2(«.7) 



-(j^^)+2(«.,j77 (j7^) 



v«^ 



+/)^> 



V«^j7^)') 



(l-2(«.4.)) 1 



(4") 



07/) 



V«^j7^)') \Jiisn"r) 

-(j7i")+2(".T37i") 07'') 



(l-2(«i4.)) 
4{(<5n,J-> 



or> 



u') 



_ (l-2(«.7) 1 
- 4{(fc,J^) r 



07/) 



V^U^ji^ V«*07/)^) 



0-+7)-2(n.tJ+7) 07/) 



- (l-2(nu)) 



/{(5j7/)^) /{(^j7/)^) 



or)-2(".tJr) 07") 



(l-2(«.4)) 1 



0^") 



07/) 



V((^j7/)^) /{(«j7/)^> 






-(j^^)+2(«,,j77 07"") 



a/((5j+/)^) /{(«j7/)^) 



or)-2(".tJi«) oi7 



/((^j7/)^) A/O^^i/T) 



(l-2(n.7) 1 
A((Sn,i,Y) r 



ij-^) 



or) 



/((^J7/)^) /((«j+/)^) 



Ini[/3, 



45j 



_ (l-2(«.7) 



or)-2(".tjr) or) 

V((^j7/)^) ^mt^r) 

0^) 



-(jr)+2("»t37") 07") 

/((«j+/)^) 4^m^) 



- (l-2(».i)) 1 

- 4((fc,^)^) r 



(^r 



/((^j7/)^) /{(«j+/)') 



and 



T^ro 1 _ (i-2(».4-)) 



-077)+2(»,,,77 0-+/) 
/((^J?/)^) AA(^7/n 



-07/)+2("'tJ7/) 0//) 



/((*J-+/)') /((^J^)'), 

or) 



- (1^2(n.;^)) 1 

- 4{(5n.^)^) r 



/((^J?/)^) /{(*J+/)') 
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